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Abstract 

The three major unknown neutrino oscillation parameters at the present juncture are the mass 
hierarchy, the octant of the mixing angle 623 and the CP phase Scp- It is well known that the 
presence of hierarchy—dc'p and octant degeneracies affects the unambiguous determination of these 
parameters. In this paper we show that a comprehensive way to study the remaining parameter 
degeneracies is in the form of a generalized hierarchy—^23 — ^CP degeneracy. This is best depicted 
as contours in the test (023 — ^cp) plane for different representative true values of parameters. 
We show that the wrong-hierarchy and/or wrong-octant solutions can be further classified into 
eight different solutions depending on whether they occur with the wrong or right value of 5 cp- 
These eight solutions are different from the original eightfold degenerate solutions and can exist, 
in principle, even if 613 is known. These multiple solutions, apart from affecting the determination 
of the true hierarchy and octant, also affect the accurate estimation of 5 cp- We identify which of 
these eight different degenerate solutions can occur in the test (023 — 5 cp) parameter space, taking 
the long-baseline experiment NOz^A running in the neutrino mode as an example. The inclusion 
of the NOi^A antineutrino run removes the wrong-octant solutions appearing with both right and 
wrong-hierarchy. Adding T2K data to this resolves the wrong hierarchy - right octant solutions 
to a large extent. The remaining wrong-hierarchy solutions can be removed by combining NOz^A 
-|- T2K with atmospheric neutrino data. We demonstrate this using ICAL@INO as the prototype 
atmospheric neutrino detector. We find that the degeneracies can be resolved at the 2 (t level by 
the combined data set, for the true parameter space considered in the study. 
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I. INTRODUCTION 


The standard three-flavor neutrino oscillation probability is described by six parameters, 
namely three mixing angles - 6 ^ 12 , 6 ^ 23 , 6 ^ 13 , two mass squared differences - Am^i 

(Amfj- = mj — m|) and the Dirac CP phase 6 cp- The neutrino oscillation data from solar, 
atmospheric, reactor and accelerator experiments have so far given information about each 
of these oscillation parameters except 6 cp [1^3]. At present, the unknowns in neutrino 
oscillation physics are (i) the sign of Amgj^ [Am‘^i > 0 known as normal hierarchy (NH) 
and < 0 known as inverted hierarchy (IH)], (ii) the octant of 623 [O 23 > 45° known as 

Higher Octant (HO) and 623 < 45° known as Lower Octant (LO)]. (iii) the CP phase Scp] 
any value of this parameter other than 0° and ±180° would signal CP violation in the lepton 
sector. In this case, it is often useful to talk in terms of the lower half-plane (LHP) with 
— 180° < Scp < 0° and upper half-plane (UHP) with 0° < Scp < 180°. 

The appearance channel P^e often known as the “golden channel” can measure all the 
three unknown parameters described above^. However, the measurement is complicated by 
the fact that different sets of values of parameters can give the same oscillation probability. 
This gives rise to degeneracies that render an unambiguous determination of true parameters 
difficult. It was discussed in Ref. [4] that there can be eight-fold degeneracies in neutrino 
oscillation probabilities which are (a) the intrinsic or 633 — Scp degeneracy [5], (b) the 
hierarchy-degeneracy [ 6 ] and (c) the intrinsic octant degeneracy [7]. The intrinsic 
degeneracy refers to clone solutions occurring due to a different 6*13 and Scp value. This 
degeneracy can be removed to a large extent by using spectral information [ 8 ]. Moreover, 
the current precision determination of 6*13 [9-12] has removed this degeneracy to a great 
extent. The hierarchy-5cp degeneracy leads to wrong-hierarchy solutions occurring for a 
different value of Scp other than the true value. The intrinsic octant degeneracy refers to 
duplicate solutions occurring for 623 and 7r/2 — 6 * 23 . 

Many papers have discussed possibilities of the resolution of these degeneracies by using 
different detectors in the same experiment [13-15]. The synergistic combination of data 
from different experiments was also discussed as an effective means of removing such degen¬ 
eracies by virtue of the fact that the oscillation probabilities offer different combinations of 

^ Originally, Rie was termed as the golden channel because of its sensitivity to 613, hierarchy and Scp- 
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parameters at varying baselines and energies [8, 16-23]. In particular, the synergy between 
long-baseline (LBL) experiments NOz/A and T2K in resolving the hierarchy-5cp degeneracy 
has been discussed recently in Refs. [24-27]. 

It has been shown in Refs. [21, 28, 29] that a precise measurement of the mixing angle 613 
is helpful for the removal of octant degeneracy. Octant sensitivity in the T2K and NOz/A 
experiments has been studied recently in Refs.[30, 31] in view of the measurement of a non 
zero 613 . The octant degeneracy is different for neutrinos and antineutrinos and hence a 
combination of these two data sets can be conducive for the removal of this degeneracy for 
most values of 6 cp [32-34]. 

Since atmospheric neutrino baselines experience strong Earth matter effects, these effec¬ 
tively remove the overlap between right and wrong-hierarchy solutions [35-38]. In particular, 
atmospheric neutrino experiments capable of distinguishing neutrinos and antineutrinos can 
be very useful in resolving degeneracies related to the mass hierarchy [38-48]. The octant 
sensitivity of the atmospheric neutrinos comes from both the appearance [49] and disappear¬ 
ance channels [50], and also benehts from signihcant matter effects, especially facilitated by 
the large value of 633 measured by reactor experiments. Atmospheric neutrinos also provide 
a synergy with LBL experiments in terms of probability behavior with respect to parame¬ 
ters, so that the combination of atmospheric neutrino data with LBL data exhibits reduced 
effect of the hierarchy and octant degeneracies [31, 45, 51-53]. 

Recently, it has been realized that for the appearance channel, the octant degeneracy 
can be generalized to the octant-( 5 cp degeneracy corresponding to any value of 6*23 in the 
opposite octant [31, 54]. A continuous generalized degeneracy in the three-dimensional 
^23 — ^13 — ^cp plane has been studied in Ref. [54]. In this work, we show that with the high 
precision measurement of ^13 by reactor experiments, the degeneracies can be discussed in 
an integrated manner in terms of a generalized hierarchy -023 - ^cp degeneracy. A good way 
to visualize the different degenerate solutions is in terms of contours in the test (023 — ^cp) 
plane for different choices of true values of parameters^. These plots also give an indication 
regarding the precision of the parameters Scp and 623 . Although hierarchy degeneracy is 
discrete, the 623 — Scp degeneracy is continuous for the appearance channel probability P^e- 

^ Note that prior to the discovery of a nonzero value of di3, the degeneracies were studied mainly in 0 i3 — 6cp 
plane. 
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Inclusion of the information from the disappearance channel restricts ^23 and discrete 
degenerate solutions are generated. We classify, for the hrst time the wrong-hierarchy and 
wrong-octant solutions with respect to right or wrong 5cp values. This also allows us to 
understand how the hierarchy and octant degeneracies can affect the precision in 5cp- We 
observe that since the wrong-hierarchy and wrong-octant solutions can occur for wrong 
values of 5cp as well, there can exist, in principle, a total of eight degenerate solutions 
corresponding to different combinations of hierarchy, octant and 5cp- This is summarized 
in Table I^. Note that these solutions are different from the eight-fold degenerate solutions 
that have been discussed in the literature. To the best of our knowledge the parameter 
degeneracies have not been studied in this generalized form in the literature prior to this. 
We identify which degenerate solutions among the eight possibilities listed in Table I exist 
in the neutrino oscillation probabilities for typical baselines and energies corresponding to 
the LBL experiments T2K and NOz/A. 

For representative true values of these parameters, we demonstrate to what extent the 
degenerate solutions can be removed by NOz/A, NOz/A-|-T2K and NOz/A-|-T2K-|-ICAL. Note 
that although the combined capability of NOz/A, T2K and ICAL in hierarchy octant and 
5cp determination have been investigated, a comprehensive study for removal of degeneracies 
using these three facilities together has not been done before. 

The paper is organized as follows. In Sec. II, we give the experimental details of the LBL 
and atmospheric neutrino experiments being considered. In Sec. Ill, hrst we summarize the 
parameter degeneracies and identify degenerate solutions at the level of neutrino oscillation 
probabilities. Then, we show their occurrence at the event level considering NOz/A and 
discuss the resolution of the different kinds of degeneracies by combinations of the given 
experiments. We also present the precision of the parameters ^23 and 5cp from the combined 
analysis with NOz/A-I-T2K-I-ICAL data. The conclusions are presented in Sec. IV. The 
Appendix A outlines the synergy between the disappearance and appearance channels and 
the role of antineutrinos. 

^ It is to be noted in this connection that if the 6 cp precision is not good then there can be continuous 
regions connecting right and wrong 6 cp solutions and hence it may not always be possible to identify 
discrete wrong 6 cp solutions. 
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Solution with 

Solution with 

right bcp 

wrong bcp 

1. RH-RO-RJcp 

V. WH-WO-WJcp 

11. RH-WO-R(5cp 

VI. RH-RO-W(5cp 

III. WH-RO-RJcp 

VII. RH-WO-WJcp 

IV. WH-WO-RJcp 

VIII. WH-RO-WJcp 


TABLE I: Various possibilities of degeneracy in the probability P^e- Here, R=right, W=wrong, 
H=hierarchy and 0 =octant. 

II. EXPERIMENTAL DETAILS 

We use the GLoBES package [55, 56] (along with the required auxiliary hies [57, 58]) 
to simulate the data of the two long-baseline nentrino oscillation experiments T2K (Tokai 
to Kamioka, Japan) and NOz/A (NnMI Off-Axis v^. Appearance, Fermilab). The sonrce 
to detector distance, L for T2K is 295 km. In the T2K experiment [59], mnon nentrinos 
are directed from J-PARC, making an off-axis angle of 2.5° toward the Snper-Kamiokande 
detector, which is a Water Cerenkov detector of mass 22.5 kt. T2K has been proposed to rnn 
for a total Inminosity of 7.8 x 10^^ protons on target (POT), and it has already collected 10% 
of the total data in the nentrino mode. At present it is rnnning in the antinentrino mode, 
and the hrst resnlts have been reported [60]. The recently operational NOz/Aexperiment is 
also sending mnon nentrinos throngh two detectors, one at Fermilab (the near detector) and 
one in northern Minnesota (the far detector), making an off-axis angle of 0.8° and traveling 
a distance of 810 km to reach the far detector, which is a 14 kt totally active scintillator 
detector (TASD). The beam power of NOz/A is planned to be 700 kW which corresponds to 
6 X 10^° POT/year which will rnn for 3(z/) -|- 3(7') years. In our simulation of NOz/A data, 
we considered the reoptimized NOz/A set np from Refs.[26, 61]. 

For the simulation of an atmospheric nentrino experiment, we consider a magnetized iron 
calorimeter detector (ICAL) planned by the India-based Nentrino Observatory (INO), the 
primary goal of which is to stndy atmospheric mnon nentrinos. ICAL@INO has an advantage 
over other detectors because it has a magnetic held which allows charge discrimination, 
thus providing the facility to distinguish between and [i~. Here, we consider the 50 


6 




kt detector for ICAL@INO with a runtime of 10 years. In our numerical analysis we have 
used constant detector angular and energy resolutions of 10° and 10% respectively and 85% 
overall efficiency. The muon resolutions from INO simulations can be found in Ref. [62]. We 
have checked that the constant resolutions used in this paper give results similar to those 
obtained using resolutions from INO simulation codes [44, 45]. In our analysis of atmospheric 
neutrinos, we used the Gaussian formula to compute the statistical . Systematic errors 
are taken into account using the method of pulls [63, 64] as outlined in Ref. [42]. We have 
added a 5% prior on sin^ 26 * 13 . 

We use the following representative values for the oscillation parameters in our numerical 
simulation as given in Table II Refs. [1-3]. 


Oscillation parameters 

True value 

Test value 

sin^ 26*13 

0.1 

0.085 - 0.115 

sin^ 6*12 

0.31 

- 

O 23 

LO=39°,42°, HO=48°,51° 

35° - 55° 

Amii 

7.60 xlO-^eV^ 

- 

Amii 

2.40 xlQ-^eV^ 

(2.15 - 2.65) xlQ-^eV^ 

Scp 

90°, 0°,-90° 

-180° to +180° 


TABLE II: Values of neutrino oscillation parameters used in our simulations. Here, the second column 
gives the true values of the parameters, and the third column represents the parameter range over which we 
have marginalized the test values. 


III. IDENTIFYING DEGENERACIES IN NEUTRINO OSCILLATION PARAM¬ 
ETERS AND THEIR RESOLUTION 

For the baselines relevant for the experiments NOz/A and T2K, the Earth matter density 
is in the range (2.3 - 3.0 g/cc), well below the matter resonance. Therefore oscillation 
probabilities computed assuming constant matter density can be used for these experiments. 
Such probabilities calculated using perturbative expansion of the small leptonic mixing angle 
6*13 (in terms of S 13 ) and the mass hierarchy parameters a (= Am^i/Arngi) are given as 
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follows [65-67], 


( 1 ) 

( 2 ) 


= 1 - sin^ 26*23 sin^ A + 0(a, Sis) 
sin^(A - 1 )A 


p — 4„2 2 

-T/ie — "±*13*23 


+ a COS 6*23 sin 26*12 


sin^ AA 


(A-1)2 

+ q;Si 3 sin 26*i2 sin 26*23 cos(A + dcp) 


A 2 

sin(A — 1) A sin AA 


(A-1) 


A 


where Sij(cij) = sin 6 *jj(cos 6 **^) for j > i {i,j = 1,2,3). We use the following notation: 
A = ArnliL/AE , A = 2EV/Amh = VL/2A, where V = y/2Gpne is the Wolfenstein matter 
term. The antineutrino oscillation probability can be obtained by replacing 6 cp ^ —6cp 
and V —V. Hence, in the neutrino oscillation probability, A is positive for NH and 
negative for IH, and for antineutrinos, the sign of A gets reversed. We observe the following 
salient features from the probability formulas: 

• The leading order term in the muon neutrino survival probability P^^, also known 
as disappearance channel, is proportional to sin^ 2023 sin^ A. This gives rise to the 
intrinsic hierarchy and octant degeneracies: 



= 

(3) 


= P/i/i(7r/2 — 023). 

(4) 


This leads to a loss of sensitivity to the hierarchy and octant, when the measurement 
is performed using this channel. 

• The appearance channel P^e does not have intrinsic degeneracies but suffers from the 
combined effect of different parameters, which leads to the following set of degeneracies: 


-P/ie(6'l3, 6cp) 

— <^CP) 

(5) 

6cp) 

= P/j.e( — ^, 6cp) ■ 

(6) 


Equations.(4, 5, 6 ) constitute the eight-fold degeneracy discussed in Ref.[4]. 

Recently, it has been discussed that in the context of probabilities which are dependent 
on sin ^023 5 the octant degeneracy can be generalized to include all values of 023 in the 
second octant [31] and can also be correlated with 6 cp [31, 68 ]. The pattern of parameter 



degeneracies in the three-dimensional ^23 — ^13 — ^cp space arising from the appearance 
probability P^e has been discussed in Ref. [54], This is a continuous degeneracy and can be 
expressed as 

PiJLe{ 021 ,-, 6^13, ^CP) = -PAte(^23) ^13) ^'cp) ^ generalized octant degeneracy. (7) 

However the reactor experiments have measured sin^ 6^13 with a high degree of accuracy and 
future measurements are expected to improve it further. This has reduced the impact of ^13 
uncertainty on octant degeneracy to a large extent [31]. In this paper, we consider another 
generalized degeneracy which is the hierarchy -023 - ^cp degeneracy: 

P^e(023, A, 5cp) = Ptie{023, “A', S'cp) ^ generalized hierarchy - 623 - Sqp degeneracy. 

( 8 ) 

This degeneracy can be observed best in the test 623 — 5cp plane. Studying it in this fashion 
allows us to view the degeneracies arising out of the remaining unknown parameters in a 
comprehensive manner. We note that while the hierarchy degeneracy is always discrete, the 
6^23 — ^cp degeneracy arising out of the appearance channel is continuous. On the other 
hand, the intrinsic octant degeneracy arising from the channel is independent of Scp 
and discrete in 623 except for 623 values close to maximal. Thus combining the survival and 
conversion probabilities gives rise to disconnected degenerate regions in the 023 — ^cp plane. 
We have elaborated on this point in the Appendix. 

In the next subsection, we study the occurrence of the above degeneracies in terms of 
probabilities for NOz/A and T2K and identify the different possible degenerate solutions at 
the probability level. 

A. Identifying degeneracies at the probability level 

Figure 1 shows the probability P^e for neutrinos (left panel) and antineutrinos (right 
panel) as a function of 5cp for both NH and IH. The plots in the upper panel correspond 
to NOz/A while those in the lower panel are for T2K. These probabilities are plotted for the 
energy where the neutrino flux peaks. The hatched area denotes variation over 023. For the 
lower (higher) octant, we vary 623 between 39° — 42° (48° — 51°). This is a good assumption 
for 023 not too close to its maximal value, for the purpose of illustrating the physics, since 
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for a given 6 * 23 (true), the disappearance channel anyway excludes values away from 6 * 23 (true) 
and 7 r /2 — 6 ^ 23 (true). Thus, these plots implicitly assume information from the disappearance 
channel. From Fig. 1, the following points can be noted. 


NOvA(E = 2.0 GeV) NOvA(E = 2.0 GeV) 




T2K(E=0.6 GeV) T2K(E=0.6 GeV) 




FIG. 1: The oscillation probability as a function of 5cp- Here, upper row represents oscillation 
probability for NOz/A [L = 810 km], and the lower row represents probability for T2K [L = 295 kmj. The 
left panel is for neutrinos, while the right panel is for antineutrinos. 

Forneutrinos— 

(i) The NH probabilities are higher than the IH probabilities. This is because of enhanced 
matter effect for neutrinos for NH in the Earth’s matter. 

(ii) For both NH and IH the probabilities are higher in the LHP. 

(iii) The probabilities for higher octant are higher for both NH and IH. 

While, for antineutrinos— 

i) “A” changes its sign, and IH probabilities become higher than NH. 
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ii) The flip in sign of Scp makes both the NH and IH probabilities higher in the UHP. 

iii) Like neutrinos, the probabilities for a higher octant remains higher for both NH and IH. 

For both neutrinos and antineutrinos the lowest line in the LO (HO) band corresponds 
to 39°(48°), while the highest point corresponds to 42° (51°), due to the sin^ ^23 dependence 
of the leading order term. 

The overlapping regions between various curves at a specihc value of the Scp indicate 
the degeneracy occurring for the right value of Scp, while the same value of probability for 
different 6 cp values denotes degeneracy occurring at wrong values of 5cp- Clearly, the former 
would correspond to solutions with the wrong hierarchy and/or octant with right 6 cp-, while 
the latter would correspond to the solutions with the wrong hierarchy and/or octant and 
wrong 5cp- The wrong-CP degenerate solutions corresponding to a given true value of 6^23 
and 6 cp can be obtained from the probability hgures by drawing a horizontal line through 
this point. The different intersection points of this line with the probability bands are 
degenerate as they share the same value of probability. However, the degenerate solutions 
occurring for a particular 6^23 (true) which is not in the vicinity of 7r/2 — 6^23 (true) in the 
opposite octant will be excluded by the disappearance channel, and the occurrence of these 
solutions in the test 6*23 — 6 cp plane will depend on the 6*23 precision of the disappearance 
channel. 

Below, we explain the occurrence of the different degenerate combinations of {hierarchy, 
O 23 , 6 cp} taking the NOz/A neutrino probabilities (the top-left panel) as reference unless 
otherwise mentioned: 

1. The overlapping regions between the NH-LO (blue) and IH-HO (green) bands around 
6 cp = — 120 ° and 90° give rise to WH-WO-R(5cp degenerate solutions in the proba¬ 
bility. However, for antineutrinos these bands are well separated. Thus, combining 
NOz/A neutrino and antineutrino data can help in removing these solutions. 

2. The probability corresponding to UHP of the NH-LO (blue) band can be the same 
as those for the LHP of the IH-LO (yellow) band for 623 = 39°. This can give rise 
to WH-RO-W(5cp solutions. Note that this degeneracy is present in the antineutrinos 
for the same values of 6 cp- Thus, for true NH, UHP (i.e. 0° < 6 cp < 180°) is 
the unfavorable half-plane of 6 cp, and this degeneracy cannot be resolved by using 
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NOz/A data alone. For T2K, the probability exhibits a sharper variation with (5cp, 
and hence this degeneracy is less pronounced between the UHP and LHP. Hence, the 
addition of T2K data to NOz/A can be helpful in removing this degeneracy. For LHP 
(i.e. —180° < 5cp < 0°), which is the favorable half-plane of 5cp in NH, there is no 
WH-RO-W(5cp solution for both NOz/A and T2K. 

3. For Scp & UHP, the NH-HO (red) band can share same value of probability with NH- 
LO (blue) band for 6 cp € LHP. This is the reason for the RH-WO-W(5cp solution. For 
antineutrinos, the degeneracy is seen to be between NH-HO-LHP and NH-LO-UHP. 
Thus, a combination of neutrinos and antineutrinos helps to remove this degeneracy. 

4. The WH-WO-W(5(7 p solution can be observed along the isoprobability line that in¬ 
tersects the NH-LO (blue) and IH-HO (green) bands at different values of 6 cp- This 
degeneracy can be seen for instance between {NH, 39°, —180°} and (IH, 51°, 0°}. 
Again, the antineutrino probability does not suffer from this degeneracy, and thus 
combining neutrino and antineutrino data can be helpful in removing these solutions. 

5. One can also have RH-RO-Wdcp solutions as a result of a so called “intrinsic CP 
degeneracy” that occurs for the same hierarchy and same value of 623 but at a different 
value of 6 cp- This is due to the harmonic dependence of the probability on 6 cp- For 
instance, within the NH-HO (blue) band, 6 cp = 0° and 6 cp ^ —135° have the same 
value of probability for 623 = 39°. However, for antineutrinos, this occurs for 6 cp = ^° 
and 6 cp = -|-135°. Thus, a combination of neutrino and antineutrino data can help to 
get rid of this degeneracy. This can also be seen to occur for T2K, for (NH, 48°, —180°} 
and (NH, 48°, 0°}. For T2K, since the flux peak coincides with the probability peak, 
the CP-dependent term is proportional to sin(5cp and thus this degeneracy occurs for 
6 cp and n — 6 cp [68]. For NOz/A, since the flux and the probability peak are not at 
the same energy, the degeneracy does not correspond exactly to 6 cp and n — 6 cp- It 
is interesting to note that this degeneracy does not occur for 6 cp = ±90°. 

Thus, among the eight solutions listed above, only the WH-RO-Rdcp and RH-WO-R(5(7p 
solutions do not exist even at the probability level. The above description is in terms of 
probabilities without including any experimental errors. At the event level, many of these 
may not appear as discrete degeneracies. In particular, for a C.L. beyond the reach of a 
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particular experiment’s precision, the different discrete degenerate solutions merge, and the 
degeneracy becomes continuous. 

Note that another way to understand the degeneracies is the biprobability ellipses in the 
P^e ~ Pfie plane [ 6 ]. This requires a single plot for neutrinos and antrineutrinos. However, 
each point on these ellipses corresponds to different values of 5cp which cannot be read off 
from the plots. The probability band plots presented in this paper provide a complementary 
way to visualize the occurrence of the degeneracies at different CP values, and one can 
readily identify the wrong and right 5cp solutions which are the main focus of this work. 

B. Identifying degeneracies at the event level 

To show the occurrence of the different degeneracies at the event level, in Fig. 2 we 
present a set of contour plots in the 6*23 — 5cp test-parameter plane assuming only neutrino 
run (6 years) of NOz/A, which is denoted as [6-1-0]. We note that the proposed run time 
of NOz/A is 3-1-3. However in this section we intend to identify which of the different 
degenerate solutions discussed in Table I can arise in the 6*23 — 5cp plane. Since the wrong- 
octant solutions disappear including the antineutrino run, the 6 - 1-0 case is the best option 
for visualizing all the possible degenerate solutions. The 3-1-3 case is discussed in the next 
section. These sets of plots also show the role of statistics which can give enhanced precision 
vis-k-vis the antineutrino run which helps in resolving degeneracies. Similar plots can also be 
drawn for the T2K 8-1-0 case. However, for T2K, since the hierarchy sensitivity is much less, 
the possibility of getting continuous regions instead of discrete degenerate solutions is more. 
Thus, the different degenerate solutions cannot be visualized so distinctly, and in this section, 
our main aim is to identify the different degenerate solutions in the ^23 plane. This can 
be done better with NOz/A 6-1-0 as the illustrative example. These plots are drawn assuming 
true hierarchy to be NH and different choices of true values of ^23 and 5cp- In this and all 
the other subsequent hgures the successive rows are for 623 = 39°, 42°, 48°, 51°. The true 5cp 
values chosen are ±90° corresponding to maximum CP violation and 0° corresponding to 
CP conservation. The blue contours correspond to the right hierarchy and magenta curves 
correspond to the wrong-hierarchy. 

The hrst column of Fig. 2 shows the degenerate solutions for 5cp = —90° for NOz/A 
running only in the neutrino mode. For 623 = 39°, apart from the true solution, RH-WO- 
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W 6 cp and WH-W0-R(5cp solutions are observed in the upper and lower right quadrants 
respectively. The RH-WO-W(5cp solution is also seen for 623 = 42°. For this case, at 
Scp = —90°, the uppermost point of the blue band in the NOz/A neutrino probability, in 
Fig. 1 one can see that the same value of probability is possible for NH-HO (red band) near 
^cp = +45° and ±180°. This explains the shape of the allowed zone - wider at these values 
and narrower at 90°. The WH-WO-R(5cp solution is seen only at a 2a level for 623 = 42°. 
This can be understood by observing that the points 42° (the upper tip of the blue band) 
and 48° (the lower tip of the green band) are more separated as compared to 39° (the lower 
tip of the blue band) and 51° (the upper tip of the green band). For 623 in the higher octant 
(48° and 51°) there are no spurious wrong-hierarchy solutions even with only neutrinos. This 
is because for NH, 623 in the higher octant and —90° correspond to the maximum probability 
for neutrinos and this cannot be matched by any other combination of parameters. Hence, 
no degenerate solutions appear and only the neutrino run for NOz/A suffices to give an 
allowed area only near the true point. Note that the contours for 48° extend to the wrong 
octant also. However, (here and elsewhere) this is not due to any degenerate behavior of the 
P^e probability but due to the poor 623 precision of the P^^ channel near maximal mixing. 

The second column represents 6cp = +90°. In this case, we observe a WH-WO-R(5cp 
solution for both 623 = 39° and 42°. This can be understood from the intersection of the 
blue and green bands in Fig. 1 close to Scp = 90° in the UHP. We also get a WH-RO-W(5 c7p 
region in the LHP. For 42°, since the 623 precision coming from the disappearance channel 
is worse, at 2 a both these solutions merge, and a discrete degenerate region is not obtained. 
For 623 = 51° from Fig. 1, we see that the point {NH, +90°, 51°} in the red band intersects 
the blue band around (NH, —90°, 39°} giving a RH-WO-W(5cp solution. 
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FIG. 2: Contour plots for NOz/A[6+0] with true values of 023 = 39°, 42°, 48°, 51° in successive rows. The 
three columns correspond to Scp = —90°,90°,0° respectively . 
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A WH-R0-W5cp solution is also obtained in this case in the LHP. Similar regions are 
also obtained for 623 = 48°. However, the RH-WO-W(5cp solution merges with the true 
solution at 2 a level. 

For 6 cp = 0 °, a discrete RH-RO-W(5cp solution is seen to be allowed at la for 623 G 
LO. This is due to the intrinsic CP degeneracy as discussed in the context of probabilities. 
But at 2a, due to the poor 6 cp precision this degeneracy becomes continuous and the whole 
LHP becomes allowed. For 623 belonging to the higher octant larger statistical errors are 
involved as compared to 6*23 G LO, and this degeneracy appears as continuous in the LHP 
even at la, and at 2a the full 6 cp range becomes allowed. For 623 = 39° and 42°, we also 
see wrong-hierarchy solutions appearing in the wrong octant. From the probability hgure, 
we identify that this degeneracy occurs around 6 cp = —30°, —180°, 180° for 623 = 39° which 
allows the LHP of 6 cp at la and the whole 6 cp range at 2a. For 623 = 42°, this degeneracy 
is seen to occur around 6 cp = —90° giving distinct degenerate solutions at the la and 2a 
levels. For 6*23 = 42°, a discrete RH-WO-W(5cp solution appears at la. From Fig. 1 , it is 
seen that {NH, 42°, 0°} has the same value of probability corresponding to {NH, 48°, 90°}. 
At 2 a, this merges with the RH-RO-W5cp solution. For 6*23 = 39°, this solution appears 
as a 2a allowed patch around {NH, 51°, 90°}. From Fig. 1, it can be seen that the above 
points are not exactly degenerate in terms of probability but due to lack of precision, they 
become allowed. For a similar reason, the 2a patch with wrong-hierarchy appears in the 
right octant for 623 = 39° and 42°. For 623 = 51° a right-hierarchy patch occurs with wrong 
octant. For 623 = 48°, because of proximity to maximal mixing, the true parameter space 
also extends to the wrong-octant. In general, we see that the CP precision is poorer for 
(5cp = 0° at this stage. This is due to the unresolved degeneracies for = 0° which lead 
to multiple allowed regions and continuous bands at 2 a. 

C. Successive resolution of degeneracies with data from different experiments 

In this section, hrst we show the status of the above degenerate regions when NOz/A runs 
in [3-1-3] conhguration. We then study the combined potential of NOz/A[3-|-3] and T2K[8-|-0] 
in resolving the degeneracies. Finally, we add the atmospheric neutrino data from ICAL 
and show its impact. 

It is well known that due to Earth matter effects and the presence of an antineutrino 
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component in the atmospheric neutrino flux, ICAL can play a prominent role in resolving 
the hierarchy and octant degeneracies [44, 62], Since there is no 5cp dependence in ICAL, 
the hierarchy and octant sensitivities are independent of 5cp- This (5cp-independent 
adds to the for NOz/A and T2K in the degenerate region, and the wrong-hierarchy and 
wrong-octant solutions can be resolved. This aids in improving both the octant and 5cp 
sensitivities of T2K and NOz/A [31, 69, 70]. Note that this is a synergistic effect and the 
combined sensitivity is better than that obtained by adding individual y^ values. 

The advantage offered by the atmospheric neutrino detector ICAL and the synergy be¬ 
tween the various experiments in removing the degeneracies can be seen from these plots. 
It is noteworthy that the allowed area in the test 6*23 — 5cp plane also gives an idea about 
the precision of these two parameters. 

Our results are presented in Figs. 3, 4 and 5 which correspond to true 5cp = —90°, 
90° and 0° respectively. In each figure the successive columns represent NOz/A[3-|-3], 
NOz/A[3-|-3]-|-T2K[8-|-0], and NOz/A[3-|-3]-|-T2K[8-|-0]-|-ICAL respectively. In these figures 
the following generic features can be noted: 

• Comparing with the NOz/A[6-|-0] panels, in all the NOz/A[3-|-3] figures we note that 
the addition of antineutrino information removes the wrong-octant degenerate regions. 
This also includes the wrong-hierarchy regions occurring with the wrong octant. For 
6^23 = 39° and 51°, the wrong-octant regions are almost completely removed. But for 
the true ^23 values 42° and 48°, both the right-hierarchy and wrong-hierarchy solutions 
extend to the wrong-octant region. 

• When T2K data is added to NOz/A[3-|-3], it helps in removing these wrong-octant 
extensions. The wrong-hierarchy right octant regions also get significantly reduced 
in size by adding T2K data to NOz/A[3-|-3]. This is due to the fact that for T2K, 
these solutions occur at different 5cp values than NOz/A. Addition of T2K data also 
improves the precision of 6^23 and 5cp ■ 

• When ICAL data are added to T2K and NOz/A, the remaining wrong-hierarchy regions 
are resolved for all the true values of 6^23 considered. The wrong-octant extensions of 
the right-hierarchy solutions are also reduced in size and the precision of ^23 improves. 
The combination of T2K-|-NOz/A-|-ICAL can resolve all the degeneracies at a 2a level 
for true 622 , = 39°, 51° for all the three 5cp values. For the 6*23 and 5cp combinations 
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of {42°, 0°} and for 48°, ±90°, 0° , there are still allowed regions in the wrong- 
octant. Note that some of the wrong-octant regions that are removed by the NOz/A 
antineutrino run could also be removed by the ICAL data. 

Apart from the above features, the following important points can be observed from the 
hgures: 

• For Scp = —90°, there are no wrong-hierarchy solutions in NOz/A[3±3], and the 
addition of T2K helps in improving the 623 precision. This feature is particularly 
prominent for 623 = 42° and 48° where T2K data help in removing the wrong-octant 
extensions for the right hierarchy solutions. With the addition of ICAL data, the 
wrong-octant solution is almost resolved for 623 = 42°, while for 623 = 48° the same is 
resolved at la. 

• For Scp = —90° and 623 = 51°, NOz/A can already resolve all the degeneracies with 
6 years of neutrino run only as can be seen in Fig. 2. However, the precision of 623 is 
worse with NOz/A[3±3]. This is because splitting the neutrino run into equal neutrino 
and antineutrino runs reduces the statistics and hence the precision becomes worse. 

• For Scp = 90°, we also see that for 623 = 48°, the wrong-hierarchy region in 
NOz/A[3±3]±T2K[8±0] is still quite large and this is where ICAL has a remarkable 
role to play. We see that when ICAL data are added, the large wrong-hierarchy region 
corresponding to 623 = 48° completely vanishes. 

• The small la wrong hierarchy-wrong octant allowed zone for Scp = 90° and 623 = 
42° in NOz/A[3±3] can be identihed as the part of the WH-WO-R5 cp solution of 
NOz/A[6±0] by comparing with Fig. 2. 

• For Scp — 0 °, adding T2K data to NOz/A[3±3] improves the precision considerably 
and also removes the wrong-hierarchy solutions to a large extent. The precision of 
6*23 and Scp around the true solution also improves. The enhanced precision due to 
adding T2K is also responsible for reducing the continuous allowed regions to discrete 
degenerate solutions for 6*23 values near 45°. Adding ICAL data removes the remaining 
wrong-hierarchy regions and further help to pinpoint the allowed zones at a 2 a level. 
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• In NOz/A[3+3], for Scp = 0° and 623 = 39°, 42°, comparing with the corresponding 
hgures in Fig. 2, we see that the spurious solution appearing at —150° at la due to 
the intrinsic degeneracy is no longer present with the addition of antineutrino data, 
since for the antineutrino probability in NOz/A the intrinsic degeneracy is between 
0° and +150° as discussed earlier. Thus the addition of neutrino and antineutrino 
data solves the intrinsic degeneracy at la at both these Scp values but at 2 a both 
±150° remain allowed. The allowed area near the true value increases in size because 
replacing half the neutrino run with antineutrinos reduces the statistics and hence the 
precision becomes worse. 

The following additional observations can be made regarding alternative parameter values 
and running modes: 

• In generating the above plots, we considered T2K running in neutrino mode with its full 
beam power. We find that once one includes the antineutrino run from NOz/A, running 
T2K in the antineutrino mode is no longer necessary for removing spurious wrong- 
octant solutions. Rather, running in the neutrino mode gives enhanced statistics and 
hence better precision. If on the other hand NOz/A runs in full neutrino mode and the 
antineutrino component comes from T2K, we have verihed that we get similar results. 

• We have presented the results for the case of true NH. If the true hierarchy is chosen 
to be IH, one would get a different set of allowed regions based on the degeneracies 
observed in Fig. 1. For example, for Scp = —90° and 623 = 39° for NOz/A[6+0] in 
the true IH case, apart from the true solution, RH-WO-W( 5 c 7 p and WH-RO-W( 5 c 7 p 
solutions would be obtained. This can be predicted from Fig. 1 (top left panel) by 
drawing a horizontal line from the bottom of the IH-LO band at Scp = —90°, which 
cuts both the IH-HO and NH-LO bands near Scp = 90°. 

The situation for NOz/A[3+3] and other combinations would be more complicated 
since the allowed regions and precision for true IH depend not only on the probability 
behavior but also on the statistics of neutrino and antineutrino data in the respective 
experiments. 

• The results are signihcantly dependent on the true value of ^ 13 , chosen here to be 
sin^ 26^13 = 0.1. Lower values of 613 (or worse ^13 precision) would lead to poorer CP 
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precision and more difficulty in removing the degeneracies. This is because 5cp is 
coupled with 6^13 in the oscillation probability. 


D. Distinguishability between 0° and 180° 

It will also be interesting to see how far the two CP conserving values 0° and 180° can be 
distinguished by the experimental setups considered. In this section, we discuss this issue. 
The true events are generated for 5cP — 0 °, 6*23 = 39° and normal hierarchy. In the test 
spectrum, we consider 5cp = 180° and marginalize over 6 * 13 . For purposes of comparison we 
also give the results for test 5cp = 90°. The results are presented in Table III. We observe 
that a 2 a sensitivity in distinguishing between 5cp = 0 ° and 5cp = 180° can be achieved by 
NOz/A+T2K. Adding ICAL data increases the sensitivity further. It is interesting to note 
that for beam based experiments, 0° and 90° have much larger separation than that between 
0° and 180°. But for ICAL, 0° and 180° are more separated though the values are very 
small. This is because ICAL itself has limited CP sensitivity due to angular smearing over 
all directions [69]. Note that for experiments like PINGU the CP sensitivity can be higher, 
and the difference between 0° and 180° can be appreciable [71]. 


1 / + p 

Test Scp = 90° 

Test Scp = 180° 

NOi/A[3+3] 

6.31 

2.82 

NOi/A[3+3]+T2K[8+0] 

14.63 

4.77 

ICAL 

1.21 

1.60 

NOi/A[3+3]+T2K[8+0]+ICAL 

14.83 

5.4 


TABLE III: sensitivity for test Scp = 90°, 180° with true 5cp = 0°, true hierarchy as NH and 

true 623 as 39°. 


E. Precision of 6*23 and 5cp 

As stated earlier, the contour plots also give information about the precision of 623 and 
5cp- In general the presence of degenerate solutions leads to a worse precision (a larger 
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width of the allowed area) in these parameters. For most values of true 6 cp and 6 * 23 , there 
is a negligible difference between the (5cp-precision of NOz/A[6+0] and NOz/A[3+3] around 
the true point. While there is a qualitative advantage to including both neutrinos and 
antineutrinos because of their different dependences on 6 cp, this advantage is squandered by 
the lower cross section of antineutrinos. The precision in these parameters can be quantihed 
using the following formulas: 

^Max _ xMin 

CP precision = x 100% (9) 

360° ^ ' 

nMax _ nMin 

023 precision = g,, g,, x 100 % ( 10 ) 

'^23 “T ^23 

In Table IV, we list the values of the la and 2a precision of 023 and Scp using these 

expressions for the case of NOz/A[3+3] + T2K[8+0] + ICAL. The CP precision is seen to be 

better for Scp = 0° as compared to Scp = ±90° for a given true value of 623 . This is because 
in the absence of degeneracies, the precision simply follows from the probability expression, 
where dP^e/dScp is smallest at 0° [25]. On the other hand, for a given value of Scp-, the 
CP precision is seen to become worse with increasing 023. The 623 precision is worse near 
maximal mixing and improves as one moves away. 


True value 

LO-Precision 

True value 

HO-Precision 

023 

Scp 

la 

2 cj 

023 

Scp 

la 

2 cj 



023 

Scp 

023 

Scp 



023 

Scp 

023 

Scp 


+90° 

1.02 

26.63 

2.17 

39.50 


+90° 

3.15 

28.45 

7.70 

48.27 

39° 

-90° 

0.89 

34.52 

2.17 

41.52 

48° 

-90° 

3.15 

30.00 

7.35 

43.22 


0 ° 

0.64 

15.83 

2.04 

28.33 


0 ° 

4.03 

17.50 

7.59 

35.80 


+90° 

1.6 

27.00 

3.32 

38.52 


+90° 

0.88 

30.32 

2.16 

43.33 

42° 

-90° 

1.7 

29.77 

3.31 

41.52 

51° 

-90° 

0.98 

34.48 

2.16 

45.00 


0 ° 

1.66 

15.83 

3.08 

29.16 


0 ° 

0.88 

19.16 

2.16 

37.50 


TABLE IV: Percentage precision of parameters 623 and Scp around true value using NOi^A+T2K+ICAL. 
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FIG. 3: Contour plots in the test 023 — 5cp plane for true 5cp = —90° and true 023 = 39°,42°, 48° and 51° 
in successive rows. The first column is for NOi^A[3+3]. The second and third columns are for NOj^A[ 3+3] 
+ T2K [8+0] and NO:/A[3+3]+T2K[8+0]+ICAL respectively. 
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FIG. 4: Same as in Fig.3 but for true Scp = +90° . 
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FIG. 5: Same as in Fig.3 but for true Scp = 0° . 
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IV. CONCLUSION 


In the era when the value of ^13 was unknown, an eight-fold degeneracy of neutrino 
oscillation parameters was identihed, which included the intrinsic ^ 13 , hierarchy-(5cp and 
octant degeneracies. With the precise measurement of ^13 the intrinsic degeneracy is largely 
removed and a four-fold degeneracy out of the original eight - involving wrong-hierarchy and 
wrong-octant solutions - remains to be solved by the current and upcoming experiments. In 
this paper, we study these degeneracies in detail and propose that the remaining degeneracies 
can be studied in the most comprehensive manner by considering the generalized hierarchy- 
^23 — ^cp degeneracy. This degeneracy is continuous for the P^e channel. The addition of 
information on the measurement of 623 by the P^^ channel gives rise to discrete solutions. 
These are best visualized by contours in the test {623 — 5cp) plane drawn for both right 
and wrong-hierarchy for different representative values of true parameters. We show that, 
depending on whether the wrong-hierarchy and/or wrong-octant solutions occur with right or 
wrong values of 5cPi there can be a total of eight possibilities. We study these possibilities 
at the probability level for T2K and NOz/A. At this level, the degeneracy is dehned as 
the equality of the probabilities for different values of parameters. However, at the 
contour level, because of the precision of the experiments, one gets hnite allowed regions 
corresponding to degenerate solutions. We dehne a degenerate solution to be one which is 
distinct from the true solution at the la level. 

Taking only the neutrino run of NOz/A as an illustrative example, we identify which 
of these degenerate solutions actually occur for different representative choices of true 
parameters. The sample true values that we consider for obtaining the contours are 
O 23 = 39°, 42°, 48° and 51° and 5cp = ±90°, 0 °. At the present level of precision, for 
Scp = ±90°, the right (wrong) Scp solutions are those which occur in the same (opposite) 
half-plane as compared to the true solution. Since 6 cp = 0° is common to both half-planes, 
for this case the right and wrong 6 cp solutions at a particular C.L. are inferred from the 
nature of the contours. The different degenerate solutions obtained are the (i)WH-WO- 
Rhcp, (ii) RH-WO-W(5cp, (iii)WH-RO-R(5cp, (iv) RH-RO-W(5 cp and (v) WH-WO-W(5 cp 
regions. Although the options i-iii have been noticed in the literature earlier, the option iv 
which exists for the same true 623 but different 6 cp has not been discussed extensively. A 
probability level discussion was done in Ref. [ 68 ], where it was called 6*23 — Scp degeneracy. 
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However, since it can occur for the same hierarchy and same ^23 we call it “intrinsic CP 
degeneracy”. The WH-W0-W(5cp solutions often appear as part of i, given the CP precision 
of the current experiments. We identify a few points in the true parameter space where this 
solution appears as a distinct degenerate solution. Note that for a true value of 623 in the 
range 48° — 51° and Scp in the lower half-plane (—180° < Scp < 0°), the NOz/A neutrino 
probability being highest cannot be matched by any other combination, and hence no de¬ 
generate solutions appear. In this case, only the neutrino run is better as it gives a better 
precision. In all other cases that we have studied, 3 years of the neutrino and 3 years of the 
antineutrino run of NOz/A are helpful in removing the wrong-octant solutions i, ii and v to 
a large extent. This also improves the CP precision since the wrong Scp solutions occurring 
with the wrong-octant are resolved. Next, we present the results combining NOz/A[3-|-3] 
with T2K[8-1-0]. It is seen that the synergy between T2K and NOz/A helps in removing the 
WH-RO-W(5cp solutions for true Scp = 0°,90°. For true 6cp = —90°, NOz/A itself is suf- 
hcient for removing this degeneracy. The precision of both parameters also improves when 
these two sets of information are compounded together. The remaining degenerate solutions 
at 2(7 can be resolved by adding ICAL data. The latter is seen to play an important role 
in removing the wrong-hierarchy solution for 623 = 48°. In conclusion, we show that the 
combination of data from different LBL and atmospheric neutrino experiments can play a 
crucial role in removing the degeneracies associated with neutrino oscillation parameters, 
thereby improving the precision of the parameters 623 and 6cp- This also paves the way 
toward an unambiguous determination of these parameters. 

Appendix A: Synergy between appearance and disappearance channel and role of 
antineutrinos 

In this Appendix, we discuss the origin of discrete degenerate regions in the test {Scp — 
623) plane from the combination of appearance and disappearance channels for NOz/A. We 
demonstrate the role of antineutrinos in resolving the degeneracies. The reference true point 
chosen in generating the data is 6cp = —90° and 623 = 39°. In the upper row of Fig. 6, 
we plot the sensitivity of NOz/A[6-|-0]. The serpentine curves in the top-left panel of Fig. 
6 represent the allowed area at 90% C.L. from only the appearance channel. The area 
inside the vertical curves represents the allowed area from only the disappearance channel 
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at the same C.L. The area between the blue dotted (magenta dotted) curves denotes the 
region obtained for the right (wrong) hierarchy. For the appearance channel, the allowed 
region is continuous and no discrete degenerate solutions appear. This can be understood 
in the following manner. In the neutrino appearance channel, 6cp = — 90 ° corresponds 
to the maximum value in the probability. As one moves away from — 90 °, the probability 
decreases and reaches its minimum value at + 90 °. On the other hand, the probability 
increases(decreases) as 623 increases(decreases). So if we draw an imaginary horizontal line 
and an imaginary vertical line at the true point, then the allowed region is expected to come 
along the diagonal of the rectangle obtained by the intersection of these two imaginary lines 
and the X,Y axes for 623 > 39 ° and 5 cp < + 90 °. For 6cp > + 90 °, the probability starts 
to increase, so 623 has to fall to keep the probability same. This explains the serpentine 
nature of the allowed area. The width of the band corresponds to the 623 precision of the 
experiment. For the disappearance channel, the allowed region is in the vicinity of 623 and 
7 r/ 2 —6*23 and parallel to the 6cp axis since the probability has a very weak dependence on 
6cp- However, the combination of the disappearance and appearance channels gives discrete 
regions in the parameter space due to the excellent 623 precision of the disappearance channel 
near 6*23 = 39 ° and 51 °. This helps to exclude the other wrong values of 6*23. This is shown in 
the top right panel of Fig. 6. Apart from the allowed regions around the true value, one can 
identify the distinct degenerate solutions corresponding to wrong hierarchy-wrong octant- 
right 6cp (WH-WO-R( 5 cp) and right hierarchy-wrong octant-wrong 6cp (RH-WO-W^cp) 
regions. 

To show the exact synergy between the appearance and disappearance channels, in the 
middle panel of the top row we plot the as a function of 6*23 (test) for a hxed 6cp value of 
— 90 ° for the same hierarchy (NH). This hgure shows that, though the disappearance channel 
suffers from the intrinsic octant degeneracy and does not have any octant sensitivity itself 
{x^ ~ 0), when added to the appearance channel, the channel is responsible for that of 6cp- 
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FIG. 6: The plots in the upper row are for NOz/A running only in the neutrino mode i.e. NOz/A[6+0]. 
Those in the lower row are for NOz^A running in equal neutrino and antineutrino mode i.e NOz/A[3+3]. 
To generate these plots, we have assumed true 023 = 39°, true 6cp = —90° and true hierarchy = NH, 
whereas test parameters are marginalized over the range given in Table II. The plots in the middle panel 
are generated for a fixed value of the test Scp (= —90°). 


Next, we discuss the role of antineutrino runs in NOuA. In the bottom row of Fig. 6, we 
plot the same figures as the top row but for the 3 year neutrino + 3 year antineutrino run. In 
the bottom left panel, we see that when antineutrino information is added to neutrino data, 
the allowed region from the appearance channel is significantly reduced. The reason is as 
follows: as 6 cp changes its sign for antineutrinos, the serpentine shape of the allowed region 
gets flipped with respect to 6 cp- This excludes the right hierarchy-wrong octant regions of 
6 cp € UHP (i.e. RH-WO-W 5 cp) and the wrong hierarchy-wrong octant regions of 6 cp G 
LHP (i.e. WH-WO-R 5 cp)- Thus, after adding the antineutrino data only the RH-RO-Rdcp 
solution remains, as can be seen from the third panel in the bottom row. From the NOz/A 
antineutrino probability figure in the top right panel of Fig. 1 , it is seen that the probability 



































of the true point cannot be matched by any points in the NH-HO or IH-HO bands. This 
means that for this true point antinentrinos are free from the degeneracies that appear with 
the wrong octant in nentrinos. Thus, the addition of antineutrino information removes the 
wrong octant solutions of NOz/A[ 6 + 0 ] that appear in the top right panels of Fig. 6. 

The nature of the disappearance channel contours are seen to remain unaltered but now 
the allowed area is slightly broader. This is because of a reduction in the overall statistics 
dne to the smaller cross sections of the antinentrinos. This is also seen in the middle panels 
where the widths of the y^contours increase. 
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